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Utilization of a quantum system whose time-development is described by the nonlinear 
Schrodinger equation in the transformation of qubits would make it possible to construct quan- 
tum algorithms which would be useful in a large class of problems. An example of such a system 
for implementing the logical NOR operation is demonstrated. 



Quantum computing algorithms 1] make use of the 
possibility of parallel operations on states which make 
up a superposed set. After appropriate operations have 
been made, a number of measurements may be carried 
out, resulting in the collapse of the wavefunction to a 
smaller set, from which information related to the solu- 
tion of the particular problem may be deduced. 

Emphasis in previous work on quantum computing al- 
gorithms has been on time independent unitary opera- 
tions on the superposed states. Focus on such operations 
is natural, as the superposed states need to be associated 
with relatively simple quantum degrees of freedom, whose 
time development could be described by a Schrodinger 
equation incorporating a simple Hamiltonian. On the 
other hand, it was noted by Abrams and LloydQ that a 
much richer set of problems can lend themselves to so- 
lution by quantum computation if nonlinear evolution of 
the qubit states could be realized. In particular, they 
demonstrate that the availability of a special nonlinear 
operation on a single qubit enables the construction of a 
quantum algorithm whose repetitive application results 
in efficient progressive separation of searched states from 
others. They also show that implementation of a two 
bit nonlinear quantum AND gate allows for an algorithm 
which finds the answers to an NP-complete problem with 
certainty in linear time. It is indicated that the two bit 
nonlinear transformation itself may be obtained through 
ordinary unitary operations in combination with single 
qubit nonlinear operations. 

Such nonlinear single qubit operations have been 
studiedH S H 0, 0, 

with the emphasis being on the 
analysis of fidelity in obtaining the results of these oper- 
ations through unitary transformations. 

A possibility for realizing nonlinear quantum opera- 
tions is through a system dynamics described by the non- 
linear Schrodinger equation which appears in the anal- 
ysis of Bose-Einstein condensation^ (BEG) and other 
contexts0. Shi[l3 has demonstrated that nonlinear 
quantum evolution is possible in Bose-Einstein conden- 
sates coupled to one another through a tunneling junc- 
tion. An explicitly nonlinear equation of motion for the 
containment well occupation coefficients has been de- 
rived. An important feature of this work is that it demon- 
strates how entanglement in BEG is realized, starting 
from fundamental considerations. The work also contains 



a discussion of the possibility of utilizing this nonlinearity 
in quantum computations. 

Admittedly, these equations involve approximations, 
and appear only because the "background" in a collective 
quantum system is treated in some mean-field form. The 
discussion of the validity of the range of the approxima- 
tion to a particular implementation may be deferred as 
a technical detail. However, the utilization of a (possibly 
macroscopic) collective quantum event as a qubit raises 
deeper questions which must be answered. We briefly 
touch upon this point in the concluding paragraph. 

In this paper, we will assume that the dynamics of 
such nonlinear Schrodinger equations may be applicable 
to operations on qubits and look into the possibilities 
introduced by such operations. It will be shown that a 
two-qubit quantum gate can be constructed with a single 
condensate (in contrast to the more than one in 10] ) in 
the presence of a non-uniform potential. It will also be 
pointed out that availability of such nonlinear quantum 
operations allow the efficient search of an optimal solu- 
tion through pairwise elimination of possibilities making 
up an extended solution set. 

A typical quantum computation algorithm utilizes the 
creation of a superposition of parallel states, usually a 
complete set of enumerable states: 

N-l 

If the qubits forming the pure state |fc > is formed 
by n quantum states with two possible eigenstates each, 
then N — 2^, and k may be taken as the number corre- 
sponding to the binary representation generated by the 
n qubits. 

The space is then enlarged to include a function of k 
in the representation: 

1^-1 >^ 1*2 >= ^ 5] |fc > |/(fc) > . (2) 

^ k=Q 

The extended space containing \f{k) > itself is made 
up of additional qubits. A measurement carried out on 
this state will yield a superposition of a smaller subset: 

|*2>-|*3>=-L=^|fc>|/(A;)> (3) 
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where the prime on the summation indicates a sum over 
a subset of the states \k > consistent with the results of 
the measurement, a total number of N' . Similar opera- 
tions involving enlargements of the space, unitary trans- 
formations on the superposed states, and measurements 
may be carried out until a final set of measurements will 
yield results relevant to the solution of the problem at 
hand. This may be in the form of a direct result yielding 
a numerical value 0|, or one may have a statistical re- 
sult, in which a sufficient number of measurements must 
be repeated to obtain an average quantity with sufficient 
accuracy 0. 

To motivate the utility of nonlinear transformations in 
quantum computing, consider the problem of searching 
through all possible ways of completing a task, to find an 
optimal one. The number k will represent one of the pos- 
sible pathways, and we will assume that the binary rep- 
resentation of this number (i.e. values of the qubits) can 
be grouped into successive "moves" or "choices" which 
must be carried out to follow this pathway. For example, 
consider the following qubit decomposition of the state 
\k >: 



Ifc >= 



gig2 y • qrn qm+l ■ ■ ■ q2m <Z2m+l ' ' ' <Z3)i 
move 1 move 2 move 3 



(4) 



This would correspond to labeling all possible path- 
ways which could be reached by a finite number of dis- 
crete moves to be chosen from a finite number of possibili- 
ties. For example, for the traveling salesman problem|l3) 
with 256 cities, one could assign consecutive m = 8 qubits 
to represent the city to be visited at that stage, and one 
would need 256 such 8-qubits to represent the complete 
trip. (Obviously, this procedure would also produce some 
pathways which are "illegal" in the way this problem is 
defined, but these will be discarded in the solution.) Al- 
ternatively, if the remaining consecutive steps in an on- 
going chess game were to be described by \k >, then 
one could use 12 consecutive qubits to describe motion 
of a piece from a general point on the 8x8 square to 
another point on the square. All possible games with 
a total of 100 moves could be represented by a total of 
12 X 100 qubits. (Again, this type of coding generates 
an overwhelming ratio of illegal moves, which need to be 
discarded.) The function \ f{k) > is arranged to hold in- 
formation about the end result of decision process, for 
example, whether the completed moves correspond to a 
"legal" sequence, and if so, what the result is. The "re- 
sult" here would be whether the game has been "won" 
or, what the total distance traveled is within the context 
of the traveling salesman problem mentioned above. 

Note that the superposed state may be factored so that 
it can be expressed as a sum over the more significant 
qubits representing the number A:, multiplying the two 



terms corresponding to the least significant qubit: 

JV/2-1 

77= E b>(|0>l/(j|0)>+|i>|/(j|i)>) 



(5) 



J=0 



where the notation (j|0) = 2j and (j|l) = 2j-|-l has been 
used. Note that for each j, there is a "preferred" choice 
between the cases (j|0) and (j|l) (based on the values of 
/(j|0) and /(i|l)), which we will label as j'. One could 
then obtain a superposed state with a reduced number 
of terms if the following transformation could be made: 



1 



(|0>|/(j|0)>+|l>|/(j|l)>)^ !/(/)> (6) 



V2 

so that the new state is 



1 



VW2 ^T^o 



Af/2-l 

E b' >!/(/) >• 



(7) 



This process could then be iterated until a single qubit 
remains, yielding its optimal value. Once this value is 
determined, the problem is reduced to the determination 
of the remaining iV— 1 qubits, for which the above process 
must be repeated, starting with these N —1 qubits. The 
number of operations necessary for the determination of 
all of the qubits then can be seen to be proportional to 
iV^, one factor of TV coming from the repetition of the 
operation for each qubit to be determined, and another 
factor from the number of transformations of the type 
shown in Eqn. 

For many problems of interest however, including the 
ones mentioned above as examples, this transformation 
cannot be achieved with unitary operations. Nonlincar- 
ity allows for "communication" between pairs of super- 
posed states in carrying out the operation in Eqn. For 
example, the simple logical NOR operation (which is re- 
lated to how the "legality" operation of the moves would 
transform) would need to have 





|0 > H 


Hi> 


|0>) - 


-> ll> 


(8) 




|0 > H 


Hi> 


ll>) - 


|0> 


(9) 




|1>H 


Hi> 


|0>) - 


|0> 


(10) 




|1>H 


Hi> 


ll>) - 


|0>. 


(11) 



The operation is clearly nonlinear, as relations 1811 II 
above demonstrate: The sum of the left sides of the ex- 
pressions |H1 and ^2 of the transformation equals the sums 
of the left hand sides of the expressions and 1101 but the 
same obviously is not true for the right hand sides. 

A similar discussion may be carried out for optimiza- 
tion problems of the traveling salesman type, obtaining 
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an optimal set of "moves" by the pairwise elimination of 
non-preferred choices. Details cannot be provided here, 
but it has to be remembered that the standard tools of 
classical logical computation (such as the logical AND 
and OR operations) may be utilized if nonlinear trans- 
formations are to be allowed. This feature then (at least 
in principle) makes accessible to quantum computation 
all classical problems which can benefit from parallelism. 

The transformation in expressions 1811 II mav be imple- 
mented as the result of time- development through a non- 
linear Hamiltonian. As an example, we will consider a 
Hamiltonian of the Gross-Pitaevskii typejQ, 0, 0| in 
which an extra potential term proportional to the square 
magnitude of the wavefunction appears. The two qubit 
state \qoqi > with qq and qi equal to or 1 will be 
taken to be related to the occupation of four sites at 
r(9o, 9i) = iloi + qij)Ax where i and j are the unit vec- 
tors in the x and y directions respectively. The system 
then corresponds to a set of four quantum sites arranged 
as a square with side Ax. The Schrodinger Equation 
describing the system will then be 



«?i^i/'('7o,9i) 



[ipiqo, qi) + ip{qo, qi) - 2ip{qo, qi)] 



2m{AxY 
\V{qo,qi) + 01-0(90, ^i)!^] VX^o, 9i) 



(12) 



where a is a measure of the nonlinearity in the system, 
and we have used the shorthand notations q = l — q and 
■0(90: 9i) = '0(r('Zo7 qi))- The first term on the right hand 
side is the kinetic energy term of the square geometry, 
it may also result through a tight-binding interpretation 
of the interaction among four quantum wells. The coef- 
ficient of this term e = ?i^/2m(Aa;)^ with energy units 
sets the physical scales of the system. The external po- 
tential V is the quantity to be "engineered" so that the 
time-development of ip has the required form. 

FigureQlshows the occupation probabilities of the four 
sites in the system as a function of time, for the initial 
conditions indicated in expressions 1811 II Note that for 
the choice of the potential values corresponding to this 
figure, the value of the ^'(1,1) component of the state 
vector yields deterministically (i.e. with magnitude ei- 
ther zero or one) the corresponding states at the right 
hand side of expressions I8I11I 

The numerical computation was performed for values 
of the unitless time parameter r = te/Ti between and 
7.665, a/e = 2.350, and the four values of the poten- 
tial F(0,0) = -0.003554e, ^(0,1) = 2.124e, F(1,0) = 
2.352e, and V{\, 1) = 0. The numerical integration of 
the Schrodinger equation was carried out by factorizing 
the kinetic and potential energy terms in the exponen- 
tial and treating each part exactly 0. This procedure 
is correct to second order in the integration time-step, 
which was taken to be Ar = 7.665 x lO""*. The four final 




FIG. 1; The squared magnitudes of -(/'(O, 0) (dashed lines), 
^(0,1) (dotted lines), 'i/)(l,0) (dot-dashed lines), and ip^l,!) 
(thick continuous lines) for four different initial conditions. 
The initial conditions correspond to the left hand terms of 
expressions 151 through fTTI from top to down respectively. Pa- 
rameters of the numerical computation are given in the text. 
Note that the final value of square magnitude of '0(1, 1) cor- 
relates with the right hand terms of expressions |H1 through 

mi 



values of V'(l, 1) are within 0.06, 0.01, 0.04, and 0.04 of 
their ideal values. Other solutions to the problem could 
be found, one needs to adjust the values of a, t, and the 
three finite values of the potential until the final value 
of "0(1,1) is within acceptable error. (The fourth value 
of the potential is the arbitrary reference of the poten- 
tial energy and was chosen as zero.) The aim at this 
stage was not to obtain a solution with overwhelming 
accuracy to an idealized model but to show that nonlin- 
ear quantum transformations are available. A physical 
realization of such a functional block would necessarily 
be more complicated and would require a more careful 
analysis. 

The implementation of the nonlinear qubit transfor- 
mation then involves the teleportation of the two initial 
qubits into the wavefunction 0, and the teleportation of 
a single qubit of information out of 0(1, 1) after a fixed 
period of time. 

In conclusion, we have shown that evolution of a quan- 
tum system with a dynamics controlled by a nonlinear 
Schrodinger equation enables nonlinear transformations 
to be carried out on qubits, which may be used to im- 
plement quantum computational algorithms with less re- 
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strictions. It may also be possible to use the nonlinearity 
to "saturate" the qubits to some ideal values close to 
their initial states, thereby implementing some error cor- 
rection. (Quantum state purification through the use of 
nonlinear transformations has been discussed in 

A question that needs to be considered at this point is 
whether the nonlinear time development is that of a true 
quantum system, as microscopic interactions always lead 
to unitary development. The assumptions that go into 
the development of the approximations that lead to the 
nonlinear Schrodinger equation may limit the applicabil- 
ity of the corresponding dynamics to the description of 
wave phenomena without the quantum features (such as 
the second sound effect in superfluid helium). It needs 
to be confirmed that the degrees of freedom described by 
this equation still maintains the indispensable quantum 
properties of interference, entanglement, and the proba- 
bility interpretation, and is not just a collective macro- 
scopic wave phenomena. Experimental work studying 
the interference effects in BEC systems 0, [20] and the 
detailed theoretical analysis of their entanglement fliij 
seem to indicate that the Gross-Pitaevskii equation may 
indeed be valid in representing genuine quantum effects 
in these systems IIJ]. Perhaps another way of looking at 
this validity is to interpret these transformations as high 
fidelity representations of nonlinearities through the use 
of the large number of quantum background degrees of 
freedom (individually obeying unitary time development) 
contained in the "mean field" of the system. 
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